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Following on from our previous work [Phys. Rev. Lett. 98, 166801 (2007)] we examine the finite 
temperature magnetothermoelectric response in the vicinity of a quantum critical point (QCP). 
We begin with general scaling considerations relevant to an arbitrary QCP, either with or without 
Lorentz invariance, and in arbitrary dimension. In view of the broad connections to high temper- 
ature superconductivity, and cold atomic gases, we focus on the quantum critical fluctuations of 
the relativistic Landau-Ginzburg theory. This paradigmatic model arises in many contexts, and 
describes the (particle-hole symmetric) superfluid-Mott insulator quantum phase transition in the 
Bose-Hubbard model. The application of a magnetic field opens up a wide range of physical observ- 
ables, and we present a detailed overview of the charge and thermal transport and thermodynamic 
response. We combine several different approaches including the epsilon expansion and associated 
Quantum Boltzmann Equation (QBE), entropy drift, and arguments based on Lorentz invariance. 
The results differ markedly from the zero field case, and we include an extended discussion of the 
finite thermal conductivity which emerges in the presence of a magnetic field. We derive an integral 
equation that governs its response and explore the crossover upon changing the magnetic field. This 
equation may be interpreted as a projection equation in the low field limit, and clearly highlights 
the important role of collision invariants (or zero modes) in the hydrodynamic regime. Using an 
epsilon expansion around three-dimensions, our analytic and numerical results interpolate between 
our previously published value and the exact limit of two-dimensional relativistic magnetohydrody- 
namics. 

PACS numbers: 73.43.Nq, 72.20.Pa, 74.25.Fy 

I. INTRODUCTION 

Quantum phase transitions (QPTs) in strongly correlated systems play an important role in modern condensed 
matter physics. In problems ranging from high temperature superconductivity to cold atomic gases, dramatic changes 
in the ground state and physical response may often accompany relatively small changes in the doping, the interaction 
strength, or other system parameters With the discovery of high temperature superconductivity in the cuprates^ 
quantum phase transitions between Mott-insulators (Mis) and supcrfluids (SFs) or superconductors, have been at the 
forefront of this intense scrutiny. More recently, remarkable advances in cold atomic gases, have allowed the observation 
of such transitions in systems of bosonic atoms^ An important stimulus for these studies, is that the SF-MI transition 
separates two of the most fascinating phases of highly correlated matter. The superfluid reveals the importance 
of phase coherence on the macroscopic scale, and the Mott-insulator the importance of strong interactions. The 
transition between the two clearly involves an interplay between strong interactions and strong quantum fluctuations. 
Such strongly correlated regimes are notoriously difficult to analyze theoretically, and shedding light on this enigmatic 
transition remains challenging. 

In this work, we focus on the magnetothermoelectric response in the vicinity of such quantum critical points (QCPs). 
Our motivation for this, and our previous work^ arose in connection with high temperature superconductivity, where 
experiments indicate strong superconducting fluctuations in a broad range of temperatures above T c £ These signatures 
appear in both thermodynamics and transport measurements performed in magnetic fields, and include enhanced 
diamagnctism^ and a strong Nernst signal^ The Nernst response is the transverse electric field induced by a thermal 
gradient in a magnetic field, and therefore hinges on the interplay of several different probes. Although a tremendous 
amount of theoretical progress has been made in various regions of the phase diagram ^ 10 ! 11 : 12 : 13 : 14 : 15 much less was 
known about the complete magnetothermoelectric response in the vicinity of such QCPsJ£ In view of the enhanced 
fluctuations, and the prospect of universal results, we advocated examining this problem at a simple, but rather 
generic SF-MI transition in the XY universality class^ — see Fig [TJ Somewhat more ambitiously, we set out with a 
view to describe the full complement of magnetothermoelectric response coefficients. 

Our strategy^ is to focus on the particle-hole symmetric transitions in the ubiquitous Bose-Hubbard model. 
This model has a distinguished history, and provides a paradigmatic example of a SF-MI quantum phase 
transitioni 17 i 18 i 19 i 20 i 21 i 22 i 23 While it arose primarily in connection with bosonic models of strongly correlated elec- 
tron systems, where the bosons are to be thought of as Cooper pairs, it has since been cleanly realized in ground 
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FIG. 1: Superfluid-Mott insulator quantum phase transition in 2+1 dimensions, as tuned by a control parameter g, such as film 
thickness or doping. The diagram shows the finite temperature 2D Kosterlitz-Thouless transition, and the crossover between 
the Mott insulator and the normal state. We examine the universal finite temperature magnetothermoelectric response in the 
vicinity of the (particle-hole symmetric) 2 + 1 XY QPT, as shown by the vertical arrow. 

breaking experiments in cold atomic Bose gases; 4 - Although not directly linked to a fermionic high temperature su- 
perconductor (with a d-wave order parameter and nodal quasiparticlcs) the simplicity of the Bose-Hubbard model 
is appealing. Amongst its many virtues, it admits a description as a quantum Landau-Ginzburg theory (or Abclian 
Higgs model) and so it naturally embraces U(l) phase fluctuations. Such fluctuations have long been argued to play an 
important role in high temperature superconductors . 24 ' 25 More recently, the finite temperature classical fluctuations 
of the Kosterlitz-Thouless^ transition have been argued to influence their diamagnetic response 14 ' 15 In this work we 
focus on the vicinity of the QCP, and study the impact of quantum critical fluctuations on the fundamental transport 
coefficients, and other physical response functions — see Fig. [TJ In contrast to other approaches, which tackle related 
problems from the superfluid side, and thus in terms of vortices, it is quite natural to examine the critical region in 
terms of particle-hole excitations of the Mott insulator. Although dual vortex formulations are possible, we do not 
pursue this complementary approach here. 

From a condensed matter perspective, our primary interests arc in the charge and thermal transport, and the 
thermodynamic response in the vicinity of the QCP. This critical point is of course well studied and there is a vast 
literature on its zero field properties which we do not attempt to review. In order to place our work in context j 5 - 
we simply recollect the most recent precursors pertaining to single field transport measurements in the absence of 
a magnetic field. It was recognized in early works that the electrical conductivity is a finite universal multiple of 
e 2 //i i 17 ' 18 i 19 ' 20 i 21 ' 22 i 27 Interactions at the fixed point are essential in order to render this finite, as opposed to a 
more conventional, non-interacting single carrier Drude peak. More specifically, this may be traced to collisions 
between counter propagating particles and holes, and in general requires a finite frequency, and finite temperature, 
hydrodynamic treatment of the problem . 20 ' 21 In contrast, it was well understood that the thermal conductivity diverges 
at the clean interacting fixed point ; 22 ' 28 ' 29 in the presence of a thermal gradient, particle and hole excitations move in 
the same direction, and the collision mechanism is unable to render this quantity finite. In the absence of a magnetic 
field, a finite thermal response therefore requires the introduction of impurities, or other forms of energy relaxation.— 

As we noted in our previous work^ this circumstance is changed markedly in the presence of a magnetic field. 
Although the DC Hall conductivity vanishes on symmetry grounds at the particlc-hole symmetric point, the application 
of a magnetic field opens up the possibility of non-vanishing thermoelectric response coefficients, even in the absence 
of impurities, or other forms of scattering. For this reason, we confined ourselves to the clean case, in order to better 
expose the main universal results. This is not a serious limitation, since results in the presence of impurities may be 
obtained from the finite frequency behavior of the clean system, provided it's not driven to a new non-trivial fixed 
point. The general problem one is interested in, is thus to apply various combinations of E, B, VT, possibly at finite 
frequency, and to measure the associated response coefficients. In view of the conceptual importance of the Bose- 
Hubbard model and the XY universality class, we present a variety of approaches to the general magnetothermoelectric 
response. 

The layout of this paper is as follows: In section [H] wc begin with a general discussion of scaling close to a QCP. 
This overview is relevant to both bosonic and fermionic systems in arbitrary dimension, and thus helps to correlate 
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the more detailed, model specific results we shall present. In section Mil we recall the field theory approach to the 
superfluid-insulator transition in the Bose-Hubbard model, and the simplifications at the particle-hole symmetric 
point. In section |IV] we describe the Boltzmann approach to quantum critical transport . 20 ' 21 In sections \V \ I VII and 
IVII1 we present a brief overview of the properties in electric fields, magnetic fields, and temperature gradients taken 
separately. In section I Villi we examine the behavior in combined crossed electric and magnetic fields, and discuss 
two distinct regimes of behavior. We approach this problem in several different ways, including entropy transport, 
Lorentz invariance, and a linear response analysis of the QBE. In section ITXl we examine the behavior in the presence 
of a temperature gradient and a magnetic field, and once again discuss two regimes. We verify that the Onsager 
relations are satisfied and obtain a non- vanishing thermal conductivity. We conclude in section 1X11 and provide several 
technical appendices. 

Whilst this longer manuscript was under construction, elegant extensions of this work appeared which also include 
the effects of impurities and particle-hole symmetry breaking at relativistic QCPs. 30 These reveal important links 
between the transport coefficients, and develop connections to the high energy community^i Applications to other 
gapless systems such as graphene have also been investigated ! 32 ' 33 ' 34 



II. SCALING FORMS 



As usual, the approach to a continuous phase transition is accompanied by a divergent correlation length, £, and a 
divergent correlation time, £ r ~ £ z , where z is the so-called dynamical exponent. Close to the transition, dependence 
on the microscopic details drops out, and non-trivial scaling relations between physical observables and the system 
parameters may be obtained. In the case of a quantum phase transition, where the microscopic energy scale is tuned 
to zero, the only relevant energy scale in the problem is the temperature, T. In this case, the divergent correlation 
time is given by 

& ~ l/T, (1) 

where we set, h = ks = lj for simplicity; sec for example the review 2 - for more details. This diverging time scale is 
accompanied by a divergent correlation length 

e~(&) 1/ *~T- 1 /* J (2) 



where a non-trivial dynamical exponent, z, reflects the potential for disparity in the temporal and static correlations. 
Throughout this work will be interested in the magnetic, electrical and thermal response in the vicinity of a QCP. 
Here, we examine the general consequences which follow from simple, yet powerful, scaling ideas. As exemplified by 
measurements on quantum Hall systems, such considerations are able to correlate a wide variety of different physical 
probes, yet are crucially independent of the microscopic approach employed^ 

The space-time dimensions of the electric and heat currents contribute to the overall temperature dependence: 

[J e ] ~ Time-Length"^ 15 ~ T 1+{d - l)/z , (3) 
[J h ] ~ Timc- 2 Length- (d - 1) ~ T 2+(d -^/ z . (4) 

In general, we are also interested in the dependence of these currents on the measuring frequency and the external fields, 
and it is temperature against which these arc compared. Having pinned the overall dimensions, these external fields 
will enter via dimensionless ratios involving the temperature^ Using the defining relations, E = A, and B = V x A, 

[E] - Length^Time -1 - T 1+1/z , (5) 
[B]~ Length" 2 ~ T 2/z , (6) 

where [A] = Length - . In this way we arrive at the following scaling forms, valid in arbitrary dimension and for 
generic dynamical exponent: 



J h (T, E, B, VT, ui) ~ T*+W* F h ( A -) , (8) 
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where, F e and F/j are universality class dependent scaling functions. One may readily incorporate additional per- 
turbations in a similar fashion. As we discuss in section IVIIi similar considerations also apply to thermodynamic 
quantities obtained from the scaling form for the free energy. We emphasize that our only assumption in deriving 
these scaling forms is proximity to a QCP. In particular, they are independent of the statistics of the underlying car- 
riers, and are equally valid for both bosonic and fermionic systems. These expressions are invaluable as they enable 
one to confirm, and sometimes infer, the field and temperature dependence of the transport coefficients. Perhaps 
more importantly, they also allow one to correlate a large number of distinct scenarios and probes, without lengthy 
or sophisticated computations. For example, as we will discuss in section [Vj in the absence of an applied magnetic 
field, or temperature gradient, linear response in E immediately yields 

J e (T, E,u) ~ T (d - 2 ^ z F(uj/T)E = <j(lu, T)E. (9) 

The general dependence of the AC conductivity on frequency and temperature is therefore easily read off i 20 ' 21 Of 
course, in order to pin the precise functional dependence on these variables, explicit calculations of the scaling 
functions are necessary, and we turn our attention to this problem below. Most crucially, at low frequencies, w C T, 
collisions at the fixed point necessitate a hydrodynamic, or quantum Boltzmann treatment of the critical regime i 20 ! 21 
For simplicity, we consider the, z = 1, relativistic field theory approach to the (particle-hole symmetric) SF-MI 
transition in the Bose-Hubbard model, but our interests, and overall approach are clearly broader. 

III. FIELD THEORY 

The Bose-Hubbard model has received considerable attention in recent yearS ) 17 ' 18 i 19 ' 20 i 21 ' 22 i 23 and describes bosons 
hopping on a lattice with amplitude t, and interacting via a short range repulsive interaction U: 

h = -tJ2( b l b J + - + ^J2 n ^ - ( 10 ) 

(ij) 1 1 

The Bose creation and annihilation operators satisfy the usual commutation relations, [6;, b^] = <5.y, where n, = b\bi, is 
the number of bosons at site i, and (i is the chemical potential. In the context of a Josephson array or superconductor, 
the bosons represent Cooper pairs of charge Q = 2e, tunnelling between superconducting regions. In general, one 
may also include the effects of disorder and long range interactions into such a model, but here we shall concentrate 
on the simplest case (fT0|) . The phase diagram of the Bose-Hubbard model is well established, and exhibits both 
superfluid and Mott insulating regions^ the latter occurring for strong enough repulsive interactions — see FigJ5] 
As a function of the chemical potential, /x, this model exhibits a series of Mott insulating "lobes" where the density 
of bosons is pinned to successive integers. At a given point within the Mott lobes, the energy cost for producing 
particle (or hole) excitations is measured by the vertical displacement to the upper (or lower) phase boundary. At 
the tips of these lobes, the energy cost to producing particle-hole excitations vanishes, and the model is particle-hole 
symmetric. In addition, the density remains constant as one enters the superfluid phase along a trajectory of constant 
chemical potential, passing through this apex. In the vicinity of these points, the SF-MI transition is described by 
the relativistic action of an interacting complex scalar field $^1 




where D = d + 1, and d is the number of spatial dimensions of the original Bose-Hubbard model (TIT))) , and the 
mass parameter, m, is set by the temperature. This is nothing but a relativistic, quantum Landau-Ginzburg action 
for the superconducting order parameter, and its associated fluctuations. It therefore represent a useful starting 
point to unravel the more general problem of phase fluctuations at SF-MI transitions. Away from these particle-hole 
symmetric points, the density changes as one enters the superfluid phase, and the action picks up an additional 
term, linear in the time derivative. Correspondingly, the dynamical exponent changes from z — 1 to z = 2; see for 
example Ref[l|. Here, we will focus on the case with z = 1, since it will allow us to employ the powerful machinery of 
relativistic quantum field theory. More general results, for arbitrary z, may be obtained by appealing to the general 
scaling arguments outlined in section ITT1 

IV. QUANTUM BOLTZMANN EQUATION 

A convenient and physically intuitive way to think about this problem, at leading order in the epsilon expansion, is 
by means of the quantum Boltzmann approach to quantum critical transport] 20 ! 21 In this framework we may regard 
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FIG. 2: Mean field phase diagram of the Bose-Hubbard model. At the tips of the Mott insulating lobes (indicated by dots) 
the system is particle-hole symmetric, and described by the relativistic quantum Landau-Ginzburg theory. 



the model (fTTj) as a gas of particle-hole excitations of the Mott insulator. The quantum Boltzmann equation (QBE) 
emerges at lowest order in the epsilon expansion and describes the impact of weak scattering (as controlled by epsilon) 
on these quasiparticlcs at the Gaussian fixed point. It takes the form of a nonlinear integro-diffcrcntial equation for 
the momentum space distribution functions, /±(k, t), of such particle and hole excitations 

^±Q(E + v k xB).^=I ± [/ +) /_], (12) 



where Vk = de^/Ok and = Vk 2 c 2 + m 2 c 4 . For simplicity we consider a spatially homogeneous system in uniform 
external fields. The collision term represents scattering between these excitations, and most crucially, incorporates 
the nonlinear interaction of the Landau-Ginzburg field theory (jlip and the associated critical fluctuations. 20 ^ 1 

T 2U ° [U ddkl ( ^-^ 

± ' 9 J t\ V 2e, 

(2^) d+1 <5(k + ki - k 2 - k 3 )S(e + ei - £ 2 - £3), (13) 
where scattering out of state k is given by 

^ ut = 2/ ± (k)/ T (k 1 )[l + /±(k 2 )][l + / T (k 3 )] 
+/±(k)/±(k 1 )[l + /±(ka)][l + /±(k 3 )], 

and we have suppressed the explicit time dependence of the distribution functions. Scattering in to state k follows by 
interchanging f± and 1 + f±. The structure of the collision term (fT5|) is readily seen by using Fermi's Golden rule, and 
represents the leading term in the e-expansion of the associated Keldysh field theory^ The factors of 1 + / remind us 
that we are dealing with a system of bosons in this example. As in standard treatments of |$| 4 theory, it is convenient 
to access the non-trivial fixed point by means of an epsilon expansion around the upper critical dimension.— At the 
Wilson-Fisher fixed point, where the renormalized mass vanishes, the bare couplings must be tuned to the values^ - 

a 4^ 2 T 2 e 247r 2 e 
- 2 = ^-, "0 = — ■ ( 14 ) 

where our spatial dimensionality is given by d = 3 — e; since we are primarily interested in two spatial dimensions, 
we shall set e = 1, at the end of any calculations. The collision term (fT5|) is therefore proportional to e 2 . As we 
shall see, the leading order epsilon expansion is illuminating both from a quantitative numerical perspective, and 
also in its ability to expose the external field and temperature dependence of physical quantities. This semiclassical 
approach is formally justified within the epsilon expansion where, at the temperatures of interest, the mean free path 
is much longer than the thermal de Broglic wavelength and the mean interparticle spacing. Further information on 
the relativistic QBE^ and applications to the quark gluon plasma, may also be found in the literature^ 
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In this Boltzmann representation, the U(l) electric current of the field theory (fTTj) takes the form 

Je = g /(S ?Vk|/+(k '^ / - (k ' t)l ' (15) 
and the heat current is given (in linear response) by 

Jh = / (^ v ^[/ + (k,i) + .MM)]. (16) 

Here, we use k to label the momentum and not wavevector; the dimensions of the currents are [J e ] = Cm-^-^s- 1 and 
[Jh] = Jm~' <I ~ 1 's" 1 . In this representation, the field theory conservation laws are related by appropriate momentum 
space integrals (or moments) of the QBE. 

The strategy is thus the same as in any application of the Boltzmann equation. We should solve the appropriate 
QBE for the non-equilibrium distribution function(s), and then compute the associated transport currents. We may 
then extract the fundamental transport coefficients defined via 



J* r \ fa a\ ( E 



J*;/ \a Rj \-VT 



(17) 



where a = Ta is required by the Onsager relations ; 38 ^ 9 ' 40 and we denote the response to a temperature gradient by 
an ovcrbar— in order to distinguish it from the thermal conductivity defined in the absence of a particle current i£i 
Knowledge of the transport coefficients is particularly useful since they allow one to correlate a wide variety of different 
measurements. For example, the Nernst coefficient, is the transverse electric field generated under open circuit 
conditions by a temperature gradient and a magnetic field. Imposing, J* r = 0, on equation (|17p one may relate this 
open circuit measurement to the more primitive transport coefficients by£ 

At the particle-hole symmetric point we have chosen to focus on, the Hall coefficient, a xy , vanishes and this may be 
reduced to 

V = \^L. (19) 

Such an approximate reduction may also occur in situations where o~ xy w 0. We see that the transverse thermoelectric 
response, a xy , is central to a discussion of the Nernst response. We shall examine this quantity in detail from several 
different perspectives, beginning in section IVIIII We shall also discuss the behavior of the other transport coefficients, 
including the thermal conductivity, R XX (B). 



V. ELECTRIC FIELD RESPONSE 



Before embarking on a detailed discussion of the general magnetothermoelectric response, it is instructive to recall 
the main known results pertaining to single field measurements ^ 7 ' 18 ! 20 As we discussed in section UH in the absence 
of any applied temperature gradients or magnetic fields, one expects the linear response electric current to behave as 

J^E,,)^^/^^^. (20) 

The corresponding conductivity therefore reads 

a(T,u;)~T^/*Z^). (21) 

For z = 1, this is borne out by the direct Boltzmann calculations of Damle and Sachdevi 2 ^ Indeed, the nontrivial 
scaling (|21[) as a function of u>/T was an important catalyst for their Boltzmann treatment; in general the limits 
uj — > and T — ► do not commute and S(0) ^ £(oo). As such the DC response at finite temperature differs from 
that obtained at zero temperature. The universal scaling function, S(w/T), may be computed numerically in the low 
frequency hydrodynamic regime by means of the QBE (|12p , and at zero frequency one obtain a 20 ' 21 



(22) 
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In d = 2, where e = 1, the DC conductivity is a universal multiple of e 2 /h: 

(2e) 2 4e 2 

ct(0) = S(0)^- ~ 1.037—. (23) 
n h 

This result is clearly appealing, and is remarkably close to both the self-dual value ; 42 ' 43 4e 2 /h, and to a number of 
early experiments on thin films.— More generally, the order of limits is also relevant to the electric field itself. A 
discussion of the non-commutativity of E,u,T, and the associated non-linear response, may be found in Refs. I22II23I . 

It is worth noting that the result (|22[) is singular in the small parameter e. This is a direct reflection of the 
proximity to the Gaussian fixed point, and that collisions must overcome the conventional Drude response for a single 
non-interacting carrier. This is possible in this two carrier system because electric current relaxation does not violate 
momentum conservation^ A steady state with a finite conductivity is therefore established. As we shall see, the 
specific leading order dependence on e is intimately tied to the physical quantity under consideration. For example, 
in section [Villi we shall demonstrate that the leading order contribution to a xy is regular in the epsilon expansion. To 
leading order, one may therefore drop the 0(e 2 ) collision term in the QBE and consider bosons of mass 0(y/e) moving 
under the action of the applied fields. This is in stark contrast to the case considered by Damle and Sachdev^ where 
the collision term was crucial in order to render the electrical conductivity finite and proportional to 1/e 2 . Nonetheless, 
as we shall discuss below, the collision term is central to a better understanding of the thermal conductivity in the 
drift regime. Here the leading order contribution begins at C(e 2 ) and is inversely related to the electrical conductivity. 
For uniformity of presentation we shall include the collision term throughout, and simplify when it is appropriate. 

VI. TEMPERATURE GRADIENT RESPONSE 

In view of the finite DC electrical conductivity discussed above, it is tempting to think that the system might also 
support a finite thermal conductivity, in the spirit of a conventional Wiedemann-Franz relation. A little reflection 
reveals that this is not the case however, and that the thermal conductivity in fact diverges in the clean system^ 
A simple way to see this is to note that under the action of the thermal gradient, particles and holes move in 
the same direction. As such, relaxation of the corresponding heat current requires energy relaxation, which is not 
contained in cither the original model (TTTj) or the Boltzmann description (fT2")h (This is in contrast to the electrical 
conductivity, where particles and holes move in opposite directions, with no net momentum. Collisions are highly 
effective in limiting the electrical current, without violating the conservation laws^i) It follows that the divergence of 
the thermal conductivity is intimately tied to the conservation of the energy momentum tensor^ in much the same 
way as happens in low-dimensional integrable systems^ As we shall discuss below, this conservation argument no 
longer applies in the presence of a magnetic ficld£ By minimally coupling the field theory (jlip to a magnetic field, 
the conservation law is modified. As such, in the presence of a magnetic field, a non-vanishing thermoelectric tensor 
and thermal conductivity may be supported, even in the clean homogeneous system^ Impurities and other scattering 
mechanisms such as umklapp processes, may of course render these quantities finite also, but in the first instance, it 
is clearly essential to understand the universal results pertaining to the clean homogeneous system. Before embarking 
on a general discussion of the interplay between a magnetic field and other probes and response functions, let us 
examine the response to a magnetic field alone. 

VII. MAGNETIC FIELD RESPONSE 

In view of the interest in fluctuation diamagnetis m 14 ' 15 it is instructive to investigate the response to a magnetic 
field alone. Recalling our scaling arguments from section [Til m the vicinity of a QCP, the only relevant energy scale 
against which to compare the magnetic field is the temperature. In the absence of electrical and thermal gradients, 
we therefore expect the free energy density to scale as£ 

^(T,B)~T 1+i /*/i(j^). (24) 

Here we have used the fact that the correlation length, £, diverges with the correlation time, £ r ~ 1/T, according to 
£ ~ (Ct) 1 / 2 , and so the prefactor is an energy density. In addition, we have used the fact that [B] ~ 1/L 2 ~ T 2 / 2 , as 
follows from our discussion in section [HI In weak magnetic fields, B <C T 2 I Z , we may expand this in powers of B: 



T(T, B) _T 1+(d - 4) / z B 2 . 



(25) 
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Here we assume that symmetry under reversal of the magnetic field ensures that only even powers of B appear in the 
weak field expansion of the scaling function fx. It follows that the linear response magnetization scales as 5 - 

M = -^-^T 1+ ^/ z B. (26) 
o B 

This is consistent with a finite temperature, diagrammatic Kubo calculation of the magnetic susceptibility of a charged 
scalar field, with m ~ T and 2=1. On the other hand, in strong magnetic fields, B 3> T 2 / z , it is natural to expand 
this in a power series in temperature, and recast the scaling relation (|24[) in the form 

T( Tl B)~B^I 2 h(j^j. (27) 

In strong fields, or low temperatures, B S> T 2 / z , we thus expect 

T{T,B) ~ B^ d+Z ^ 2 . (28) 

In the relativistic problem, this yields the strong field behavior, T ~ _B 3 / 2 , in d = 2, and T ~ B 2 , in d = 3. Broadly 
speaking, this non-trivial field dependence is a reflection of the relativistic Landau level spectrum, where £ ~ \[B. 
This is suggested on dimensional grounds by equation (|24[) for 2 = 1, and is borne out in a direct computation of the 
partition function of a free massive relativistic charged scalar field i 45 ' 46 These considerations are also compatible with 
elegant and highly non-trivial results obtained in the early days of quantum electrodynamics on vacuum polarization 
and pair productionj 47 i 48 i 49 see for example Ref. HO for a recent review of Heisenberg-Euler effective Lagrangians in 
both spinor and scalar quantum electrodynamics^ These detailed studies require ultraviolet rcgularization and yield 
an additional logarithm at strong magnetic fields, so that T ~ B 2 h\{l/B) in d = 3. Although this is not captured 
by the simple scaling arguments, it is intimately connected to the vacuum screening properties of the field theory via 
the renormalization group beta function. 

In this section we have discussed the distinct field regimes of the static magnetic response. In the subsequent 
discussion, we shall see how such regime divisions also emerge in the transport properties. In many respects this is 
rather natural, since the dynamics and statics are intimately related at QCPsi^ 



VIII. CROSSED ELECTRIC AND MAGNETIC FIELDS 



In order to understand the quantum critical transport in combined E and B fields, it is instructive to recall the 
motion of a single relativistic particle in crossed electric and magnetic fields^ This will provide significant orientation 
for the more general interacting field theory. 5 The most important feature is that the motion of a single charged particle 
is qualitatively different depending on whether, E < cqB, or E > cqB, where cq is the speed of light. This is most 
easily understood from the vantage point of a moving frame of referenced In the former case, there always exists a 
moving frame where the electric field vanishes, and the particle experiences a pure magnetic field. Conversely, in the 
latter case, there always exists a frame where the magnetic field vanishes, and the particle experiences a pure electric 
field. Boosting back to the lab frame, we mix in the complementary field component, but the qualitative character 
of the motion is more "electric field like" or "magnetic field like" , as dictated by the inequalities — see Fig. [3] We 
shall discuss this in more detail below, but before doing so, we note that our problem is a little more subtle due 
to the appearance of an effective speed of light, c, in the effective field theory ifTTj) , and the associated QBE p2[) . 
What matters for our purposes, at least in the first instance, is the solutions to the QBE (fT2"|) for a given fixed ratio 
of E and B. It may be seen that these solutions are kinematically distinct for, E < cB, and E > cB, as follows 
directly from the left hand side of the differential equation, without recourse to Lorentz invariance arguments. As, 
such it the effective Lorentz structure which plays a key role in determining the solutions of the relativistic QBE (fT2"j) . 
and we shall henceforth use this effective speed of light in our subsequent discussions. It is interesting to note that 
such Lorentz transformation arguments also find applications in other effective relativistic systems, as was recently 
discussed by Lukose et al in the context of graphene.— 

Under a Lorentz transformation with a velocity v, the electric and magnetic fields transform according to 5 -^ 

E' = 7 (E + v x B) - 2— ^ L , (29) 
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where, 7=1/ y/l — v 2 /c 2 , and we work in SI units. It is readily seen that for crossed E and B fields, with |E| < c|B|, 
there exists a frame moving at the drift velocity 

ExB 

v ° = iBp"' ( } 

where the electric field E' vanishes. In this moving frame, the particle is subject to a pure magnetic field of reduced 
strength, B' = B/7D. Taking our magnetic field to point along the z-axis, the particle executes cyclotron orbits in the 
moving frame with, x'(t') = rcos(u>'t') and y'(t') = =Fr sin(u/i'); here r is the radius of the orbit, 0/ is the cyclotron 
frequency, and the signs indicate the sense of rotation for positive and negative charges respectively. With the electric 
field along the x-axis, the drift velocity points along the negative y-axis. Boosting back to the lab frame using the 
inverse Lorentz transformations one finds the parametric equations of motion: 

x{t) = rcos(wY), y(t) =y(Trsm(u)'t') -v D t'), t = j(t' ± v B r sin{uj' t')/c 2 ). (32) 

In the nonrelativistic limit c — > 00, these reduce to the parametric equations of a trochoid: x 2 (t) + (y(t) + vut) 2 = r 2 — 
see Fig [3] In view of the cyclotron motion of equal numbers of particles and holes, it follows that the DC conductivity, 
&xx{B) = 0, in this regime, at least at the single particle level. On the other hand we see that the thermoelectric 
tensor, a xy , may be finite due to the finiteness of the drift velocity. We shall see that these expectations are borne out, 
even in the presence of interactions at the fixed point, and that this single particle description captures the relevant 
physics^ It also underpins the divergence of the Nernst coefficient (fT9")l in the clean, particle-hole symmetric case£ 

Returning to our single particle problem, for crossed E and B fields with |E| > c|B|, there exists a frame moving 
at velocity 

,/ExB\ 

vb = c 2 ^ J , (33) 

where the magnetic field B' vanishes. In this moving frame, the charged particle is subject to a pure electric field 
of reduced strength, E' = E/7B. In the absence of scattering, the energy e', and the components of momentum p' 
parallel to E, continue to increase indefinitely. In the lab frame 

e = 7B(e' + VB.p'), (34) 
P|| = 7B(pj| +vb£'/c 2 ), (35) 
P± = P±, (36) 
where the labels parallel and perpendicular are with respect to the boost velocity, vb . In these notations p{, is constant 



(since it is transverse to the electric field) whilst p' ± and s' = \/ c 2 p' 2 + m 2 c 4 increase with time. It follows from 
equations (|34[) , (|35p and (|36[) that the energy, and both components of the momentum increase indefinitely in the lab 
frame. That is to say, if we apply an electric field |E| > c|B|, both particles and holes acquire an identical and ever 
increasing component of the momentum at right angles to the electric fields This cannot be relaxed by the collision 
term. Since, e^k = c 2 k, we expect that a xy is infinite in this regime. On the other hand, compatibility with the 
results of Damle and Sachdcv, for B = 0, suggests the possibility of a finite value of a xx . Once again, the underlying 
divergence of the Nernst coefficient (fl~9|) is apparent. 

From the preceeding discussion, we see how the ratio, E/B, may influence transport measurements. As we shall 
discuss in section HXl there are analogous regimes in a thermal gradient, for VT < B, and VT > B, at least within 
the framework of the collisionless Boltzmann equation with a linearized driving term. Although we no longer have 
Lorentz invariancc arguments, the distinction once again shows up in the single particle kinematics of the associated 
QBE, since the thermal gradient acts like a momentum-dependent electric field. In the remainder of this section, 
we shall take the electric field regimes in turn, and examine the magnctothermoelectric transport coefficients from a 
variety of different perspectives. We will focus primarily on the drift regime, |E| < c|B|, since it is both the simplest 
to analyze, and also pertains to conventional linear response measurements at fixed magnetic field. 



A. Drift Regime: |E| < c|B| 

In this section we are interested in the thermoelectric response in the drift regime. It turns out that there are 
several complementary ways to address this problem, and we examine these below. Our strategy is to first develop 
an understanding of the QBE (fT2"|) as written, and to defer discussion of the epsilon expansion itself until later. This 
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FIG. 3: Motion of a single relativistic charged particle in crossed electric and magnetic fields, (a) In the regime |E| < c|B| the 
particle executes cyclotron orbits and has a well defined transverse drift velocity vd- 
is continually accelerated by the electric field. 



(b) In the regime |E| > c|B| the particle 



combination of different perspectives is particularly useful in establishing the Onsager relations between the trans- 
port coefficientSi 38 ' 39 i 40 It also helps demonstrate the equivalence between the field theory and quantum Boltzmann 
approaches. We begin in section fVIII A II with a computation of the thermoelectric tensor, a xy , based on the entropy 
drift of a charged scalar field. In section Will A 21 we turn our attention to the more general problem of the Lorentz 
invariant solutions of the QBE. In section IVIII A 3l we use this explicit distribution to determine the heat current 
response to an electric field, and verify the answer obtained by entropy drift. In section fVIII A 41 we demonstrate how 
this result also follows from a more familiar linearization of the QBE. Such a linearization approach will be particularly 
useful in the presence of thermal gradients, where Lorentz field transformation arguments are not available. In section 
IVIII A 5l we finally turn to the epsilon expansion itself, and evaluate our general expression for a xy numerically. We 
compare the results to our scaling analysis of section [Til We provide a very brief discussion of the non-drift regime in 
section [VaTS 



1. Entropy Flow 

As follows from the defining relations (fT7|) . the thermoelectric tensor, a, quantifies the electric current which flows 
in response to a temperature gradient. Provided the Onsager relations are satisfied, this may also be obtained (up 
to a factor of temperature) from the heat current which flows in response to a temperature gradient. The latter 
route is easier to begin with, and we consider the complementary approach in section HXl In view of the well defined 
drift velocity it is natural to compute the transverse thermoelectric response, a xyi as a transport of entropy with the 
characteristic flow 

E x B 

vd = -j^p (37) 

To lowest order in the epsilon expansion, it is sufficient to compute the entropy density of a free massive charged 
scalar field, where the mass parameter is given by (|14p . That is to say, the mass incorporates the leading Hartree 
contribution of the self-interactions ^ A free boson is essentially a harmonic oscillator, and the partition function of 
the latter is readily seen to be 

Z = Y e -/3(n+V2)^ = 1 (38) 

h 2sinh(^) 

Taking the logarithm and integrating over all momenta, the analogous expression for a massive neutral scalar field 
follows immediately^ 3 - 



where fwJk = £k — \fk 2 c 2 + m 2 c 4 . Dropping the first (divergent zero point) contribution, the free energy density reads 

T=kBT S^k in{i - e ~^- (4o) 



11 



The corresponding entropy density is given by 

(3e k 



dT , f d d k 
<i> = — -n — = —kn 
dT 



J {^h) d [ v > 



(41) 



Performing an integration by parts on the logarithmic term one obtains 

1 f d d k 



S dT 



— /o (e fe )V k .(e fe k), (42) 



where /o = (e^ efc — 1) 1 is the Bosc distribution function, and we have used the identity V k -(£fck) = v k .k + de k . 
Since particles and holes acquire the same drift velocity, the corresponding heat current density for our charged scalar 
field theory is given by Jh = 2T5vd- In particular, if we apply an electric field, E x , in the positive rr-dircction, and 
a magnetic field, B z = B, in the positive z-direction, the drift velocity points in the negative y-direction: 

jv = -2TS^. (43) 



Assuming the validity of the Onsager relations, which we demonstrate are satisfied in section UX\ we thus obtain 

a yx = --g. (44) 

Equivalently, using the Onsager symmetry relation^ 0^(13) = a yx {— B) 

25 2 f d d k 



2 /o(£fc)V k .(e fc k). (45) 



u B dBT J (2nh) 

A useful alternative form of this result, which better exposes the relation to heat currents, is obtained by performing 
another integration by parts: 



\ £fc k.(-V k io). (46) 



Using the identity 



we may also write this in scalar form as 



dBT J {2nh) 



2S 2c 2 f d d k f df 



B dBT J (2irh) d V de k 

Note that the result for a xy is positive in sign; an electric field in the y-direction and a magnetic field in ^-direction 
yields E x B drift, and thus a transport of heat, in the positive a;-direction. In the next sections we shall see how this 
result (|48|) emerges from the QBE (|12[) . even in the presence of interactions. Developing the Boltzmann approach is 
particularly useful since it will allow an intuitive and systematic calculation of all the transport coefficients in this 
regime. Before closing the section let us note that the entropy density per species of carrier has a particularly simple 
form in the massless limit. With m = one obtains 

S=C d {k B X T d ) , (49) 

where At = frc/ksT is a thermal wavelength for massless particles, and 

I [ d d k Pe* 

Ld ~ dj (2^^-1)2' ^ 

where k = ck/ksT are dimcnsionless variables. In particular, C3 = 27r 2 /45 ~ 0.439, and C2 = 3£(3)/27r w 0.574. We 
shall employ the first of these results within the cpsilon expansion in section IVIII A 51 
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2. Lorentz Invariance and the QBE 

The first step in developing the Boltzmann approach is to solve the QBE for the non-equilibrium distribution 
functions, f±(k,t). In the absence of any electric and magnetic fields, it is readily verified that the QBE (TT2")) . is 
satisfied by the static Bose distribution function 

f ± (k,t) = f (s k ) = - j ^— J , (51) 

even in the presence of the collision term. This is easily seen by noting that 1 + /o(£k) = e ■ /o(^k) ■ It follows that 
equals F± provided e(k) + e(ki) = e(k2) + e(ka). The energy conserving delta function present in equation (fT3|) 
is therefore enough to ensure that the collision term vanishes and that the full QBE is satisfied. 

Let us now examine the non-equilibrium situation in the presence of crossed electric and magnetic fields. As 
discussed in section l"VIII[ in the regime |E| < c|B|, we may move to a frame with velocity 

E x B 

VD = IbF' *■ ' 

where the electric field vanishes. Since a pure magnetic field does not affect the energy of a charged particle, it follows 
than an equilibrium distribution /o(£ k ) must satisfy the Boltzmann equation in this boosted frame. (Indeed, it is 
readily verified by direct substitution that /±(k) = /o(e k ) is a solution of the original QBE (|12[) when E = and 
B^O, since dfo/dk = v^d/o/ctek-) Since 

£ k = 7D(£k - v D .k), (53) 

where, 70 = (1 — v^/c 2 ) -1 / 2 , we conclude that 

/±(k) = / (4) = fo ( 1=A ) > (54) 

W 1 - v d/ c / 

is a solution of the full QBE (fP2|) . Again, it may be verified by direct substitution that (|54|) is in fact a solution of 
(fT2|) in the presence of our crossed E and B fields. Explicitly, the left hand side of the Boltzmann equation is readily 
seen to give zero using the fact that B x vd = E: 

(E + v k x B).^± = (E + v k x B).(v k - v D ) 7 D^4p^ oc v k .(E - B x v D ) = 0. (55) 

OK OE k 

Likewise, in the collision term we see that F±- t equals f± provided that e'(k) + e'(ki) = e'(k 2 ) + e'(k 3 ). Since 
e'(k) = 7d(£ic — VD-k), this is ensured by the combined action of the energy and momentum conserving delta 
functions appearing in (|13|) . In this way we have established the non equilibrium solution (|54p of the QBE valid in 
the drift regime |E| < c|B|. This is analogous to Kohn's theorem^ for non-relativistic electron systems.— In the next 
section we shall use this distribution to compute the transport coefficients. 

3. Transverse heat current in response to an electric field 



The equality of the particle and hole distribution functions (|54|) reflects the fact that drift velocity (|52|) is independent 
of the charge of the carrier. It follows immediately from the definition (|15[) . that the both the longitudinal and the 
transverse components of the electrical conductivity vanish in this limit. Note that there is no conflict with the 
existing results of Damle and Sachdev , 20 ' 21 since our present results are derived in the drift regime |E| < c|B|. As 
such we cannot simply set B = and recover the results obtained in the absence of B. Moreover, the vanishing of 
the conductivities are consistent with the single particle picture presented in Fig. [3] In contrast, it is clear from the 
definition (|16ll . that a non- vanishing heat current may be supported in crossed E and B fields. Substituting the drift 
solution ([54]) into (fT6]) yields 

J» = 2 / 7^7 c 2 k/ ( £k ' VD ' k ) . (56) 
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Taylor expanding the distribution function in powers of vd gives 

d d k 



{2TTH) d 



/oM+ (e / £k _ 1)2 k-v D +0(vg ) ) 



(57) 



where the first term vanishes upon integration. In the presence of an electric field, E x , and a magnetic field, B z = B, 
the drift velocity vd is in the negative y-dircction. As such, to linear order in the electric field 



T y - -2 



It follows from the defining relations (fT7|) that 



d d k 
(2nh) d 



d d k 



2;„2 



2 7 2 

■ c k 



V ( e Pe k _ !)2 B 



1 



Equivalently, using the Onsager relations^ a yx (B) = Ta yx (H) and a xy (B) = a yx {— B) 



2c 2 



d d k 



( e fSe k _ !)2 Bz 



dBT J {2Trh) d 



k' - 



del 



(58) 



(59) 



(60) 



This coincides with our previous result ([48]) obtained by entropy drift arguments. This supports the validity of our 
simple result, that at the strongly fluctuating SF-MI transition, a xy = 2S/B, where S is the total entropy density 
pertaining to each carrier type£ 



4- Linear Response 

In order to go beyond our existing results, especially with a view to thermal gradients in section IIX1 it is useful 
to perform a systematic linear analysis of the QBE. In the presence of electric and magnetic fields the Boltzmann 
equation reads 

|iQ(E|v k xB).| = I ± [/ + , /_]. (61) 

To linear order in the electric field we may parameterize the departure from equilibrium 

/±(k) = f (e k ) ± Qk.E^(fc) + k.(E x B) (62) 

where for simplicity we focus on static solutions. Here, ip(k), represents the longitudinal shift due to the applied 
electric field, and ip±(k) represents the transverse shift due to the combined electric and magnetic field. Substituting 
this expansion into (|61[) and retaining only terms up 0(E) 



±Qv k . 



( H ) 



L±[±Qk.EV] +L' ± [k.(E x B)V_l], (63) 



where the terms on the right hand side correspond to the distinct linearizations of the collision term arising from 
the longitudinal and transverse shifts of the distribution function — see Appendix [A] These are straightforward, 
but rather tedious to derive, and involve momentum space integrals over the remaining products of Bosc distribution 
functions. As noted by Damle and Sachdev ) 20 i 21 these integrals may be evaluated in terms of polylogarithm functions. 
Since the explicit form of these linearizations was not included in any of the original publications^ ' 20 ! 21 we provide 
them in Appendix [A] Although the details of the collision term are certainly important for a quantitative numerical 
implementation of the epsilon expansion, the physical results are primarily determined by the robust symmetry 
properties of these linearizations. In the case at hand, to lowest order in the epsilon expansion the collision terms 
may be neglected in accordance with our previous findings and we obtain 

Mk)=0, Mk) = ^io(~). (64) 



B| 2 V ds k 

The vanishing of the longitudinal shift is consistent with the vanishing of the DC electrical conductivity, a xx (B), in 
the drift regime. The heat current is readily found from the expansion (|62[) 



^ (2^ £fcVfek - (ExB) ^ (A:) - (65) 
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As such, to lowest order in the epsilon expansion 

j ^/(^ k < k ->(-t)- (66) 

It is readily seen that this is consistent with our original transport coefficient (|48[) . Indeed, our lowest order linear 
response solution (|64|) . yields 

/±(k) = /o(£ k )-k.v D (jg). (67) 
This coincides with the linearization of our exact boost distribution (|54[) . 



5. Numerical Evaluation of a xy 



Having presented compelling evidence for the general form of the thermoelectric tensor, let us finally obtain the 
explicit numerical result within the epsilon expansion. Performing the angular integrals yields 



t-xy 



dBT{2Trh) d J 



dk k 



d+l 



j3e 



0S1, 



( e /3e k _ !)2 ' 



(68) 



where S d = 2n d / 2 /T(d/2) is the surface area of a unit hypersphere in rf-dimensions. Introducing rescaled variables 
k = ck/k B T and m = me 2 jk B T we find 



k B ( k B T 



B V he 



2S d 



d(2ir) d J 



dk k 



d+l 



(1 



+m 2 \2 



(69) 



Within the epsilon expansion, d = 3 — e, and the mass parameter entering equation ([69"]) is proportional to ^/e as 
indicated in p4[) . To leading order we may thus evaluate the dimensionless numerical prefactor in three dimensions 
so as to obtain 



2ki 
~~B 



b I 



knT 



he 



3-e 



1 



4„-fe 



dk 



fc 4 e 



6tt 2 J (1 - e" fc ) 2 



The term in square brackets is just our massless entropy prefactor, C3 = 2ir 2 /45, and we thus find 



4tt 2 k B (k B T 



•*xy 



In particular, in d = 2, one obtains 



45 B he 



T 2 

a xy = a [ — 



3-e 



(70) 



(71) 



(72) 



where ao ~ 0.88 k B (hc)~ 2 . It is readily seen that both answers (|71|) and (|72|) are consistent with the general scaling 
arguments presented in section [Til Noting the inverse magnetic field dependence inherited from the drift velocity (|31[) . 
it follows from equation ([8]) that the heat current scales as 



j ._. T 2+(d-l)/z g 1 



T 



rpd/2 



E = &E, 



rpx+i/z b 

in response to applied electric and magnetic fields. Our results arc therefore consistent with the scaling behavior 

T d/ 



(73) 



<-xy 



B ' 



(74) 



when specialized to the Lorentz invariant case with z = 1. This completes our initial survey of the drift regime in 
crossed electric and magnetic fields. 
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B. Collision Dominated Regime: |E| > c|B| 



Having discussed the drift regime in considerable detail, let us very briefly comment on our expectations when 
E > cB. As we discussed in section [Villi a single particle is continually accelerated by the electric field. Since the 
collisions in the QBE \12\i conserve energy and momentum, there is no way to relax a transverse heat current, and so 
we expect a xy to diverge at the clean fixed point in the absence of impurities. On the other hand, the regime E > cB 
is continuously connected to that studied by Damle and Sachdev ) 20 ' 21 with E ^ and B = 0. We thus expect the 
possibility of a finite electrical conductivity, a xx (B). It would be interesting to check these expectations by including 
impurity scattering in the framework of the QBE, although we do not pursue this here. 



IX. TEMPERATURE GRADIENT AND MAGNETIC FIELD 



Having discussed the heat current which flows in response to an electric field, we now turn our attention to the 
complementary problem of the electric current which flows in response to a temperature gradient. This is important 
in order to verify the Onsager reciprocity relations , 38 i 39 i 40 which we have used at several points to recast a xy = Ta xy . 
Although we no longer have the luxury of Lorentz invariancc arguments, we will again argue in favor if two distinct 
regimes of behavior. We will begin in section IPC Al with a discussion of the drift regime. In section IIX A II we will 
recover our previous results for a xy , by taking an appropriate moment of the QBE in a thermal gradient. In section 
IIX A 21 we will further rederive this from a linearization of the collision term. Both of these approaches indicate the 
presence of a finite thermal conductivity, R XX (B), and we will examine this in more detail in section I"IX A 31 



A. Drift Regime: VT < B 

1. Transverse electrical current in response to a temperature gradient 

Although we have obtained a xy by several different routes in crossed electric and magnetic fields, it is prudent 
to consider the computation in a thermal gradient. The celebrated Onsager relation a 38 i 39 ' 40 tell us that we ought 
to be able to compute a xy by looking at the transverse electrical current which flows in response to a temperature 
gradient. Although this Onsager relation is expected to be true, it ought to be verified by direct calculation. An 
additional motivation to examine this complementary approach is that it will also pave the way to an analysis of the 
thermal conductivity, R(B). Let us therefore drop the electric field from the Boltzmann equation (|12|) and introduce 
a temperature gradient. This is easily achieved by allowing the temperature variable to be a function of position^ 
The generic form of the Boltzmann equation with E = reads 

^+v k .^ ± Q(v k xB).^=I ± [/ +J /_]. (75) 

In the absence of any material inhomogcneity we may assume that any spatial variation is due to the imposed 
temperature gradient: 

Within linear response in VT, we may replace f± by /o in the spatial gradient term. In this way we may write 

+*Vj(-£g)±«<^B).^- W+ ./_,. (77) 

A characteristic feature of the (single particle) drift regime is that there is no net energy gain during each revolution; 
the energy gained as a particle is accelerated under the electric field is lost on the reverse journey^i This suggests that 
we ought to look at the time variation of the total energy current within Boltzmann theory. We may do so directly 
from equation (|7T|) by multiplying both sides by £k"Vk and integrating over all momenta. For relativistic particles this 
is aided by the fact that e^v^ = c 2 k is proportional to the momentum: 
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where we define U = (— VT)/T, and are the energy currents pertaining to particles and holes. In writing this 
equation we have used the fact that the integral involving the collision term vanishes; the momentum is an example of 
a so-called summational invariant.— These quantities play an important role in the hydrodynamic limit. Integrating 
the third term by parts one obtains 



where we have used the fact that the terms involving derivatives of vanish. Adding the particle and hole solutions 
yields 

^ + 2c 4 /A k(kU) ^ 
dt +2C J (2nh)d k[k - [J >[de l 



2c / To — few k ( k - u ) hr^ - c J - x B = 0, (80) 



where J e and J e are the total heat and electrical currents defined by equations (|15p and (|16p . If we impose the time 
independence of J^; we obtain 

J,xB = 2c*/^( k .U)(g). (81) 

In particular, if we apply a temperature gradient, U y , in the y-direction and a magnetic field B z = B in the z-dircction, 
one obtains a transverse electrical current 

2 1 d d k , 2 ( d fo 



J * B=2c iww* K \-& ]Uy - (82) 



From the defining relations ()17|) we thus obtain 

2c 2 f d d k 2 / df 



'■HI 



dBT J {2irh) d V d£ k 



(83) 



This is in agreement with our drift answer (|48[) and our complementary calculation (|60[) based on Lorentz invariancc 
arguments in crossed E and B fields. By taking appropriate moments of the Boltzmann equation^ we have obtained 
the thermoelectric tensor without recourse to an explicit solution. Moreover, we have recovered the correct Onsager 
symmetry relation. In addition, the vanishing of J^; indicates that the thermal conductivity is finite in this regime. 
We shall examine this further in section HX A 31 Before doing this we first need to perform a linear response analysis 
in VT. 



2. Explicit Construction of Linear Response in VT 
In the presence of a temperature gradient and a magnetic field 

where U = (— VT)/T. To linear order in the temperature gradient we parameterize 

/±(k) = /ofe) + k.U#)±Qk.(UxB)^(fc). (85) 

Note that ^(fc) and tp±(k) are distinct from those defined previously in equation (|62|) . Substituting (f85|) into (|84|) and 
retaining only those terms of C(U) 

= L' ± [k.UV>] +L±[±Qk.(U x B)V_l], (86) 



c 2 k. 



U|^U(UxB)#) + «^) 



where L and L' are the distinct linearizations of the collision term discussed in Appendix [X] To lowest order in the 
cpsilon expansion we may drop the collision terms to obtain 
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The electrical current obtained from (|85|) yields 

J e = 2Q 2 J ^L Vfc k.(U x B)^x(fc). (88) 

Substituting ((HZ) into {88} we find 

7 (27rft) d V |B| 2 J \ de k J K ' 

This has a structure that is close to that of equation (|66[) . and once again this yields our previous expression for a xy . 
However, in addition we will be able to go beyond this result and examine the thermal conductivity, R XX (B). 



3. Thermal Conductivity 

Thus far, we have used the QBE in a thermal gradient and a magnetic field to verify our previous result for a xy , 
and the Onsager relation. As we have seen, this is encoded in the lowest order 0(e ) solution given in equation (|87[) . 
However, self consistency of this non-trivial result for -0j_(fc), in turn induces a non-trivial longitudinal shift (and a 
finite thermal conductivity) at the higher order of 0(e 2 ). As indicated by equation (|86|) . the O(e ) solution for ip±(k) 
yields a non-trivial tp(k) at 0(e 2 ): 

TQv t .(U x B)V(fc) = L±[±Qk.(U x B) tp±(k)] . (90) 

The right hand side of this equation is structurally similar to the linearization encountered by Damle and Sachdev^ 
in the context of (zero magnetic field) electrical transport. The only crucial difference is that the electric field has 
been replaced by U x B, where we recall that U = (— VT)/T. Using their notation, equation ([90]) may be written 



Tv fc .C^(fc) = ±k.c|-e 2 (^)^ 



dki [V>±(fc) F!(fc, h) + V>±(fci) F 2 (fc, h)} \ , (91) 



where we denote the combination Q(XJ x B) = C — see equation (|A9|) . Here, Fi(fc, ki) and F2(fc, ki) are non-trivial 
kernels which we provide in equations (|A16|) . (| A17|) and (|A18|) of Appendix [5J In writing equation (|91[) we have also 
restored a factor of c/h which stems from requirement that the collision term have dimensions of s , the kernels are 
dimcnsionlcss and the prcfactor combines with the measure of integration. That is to say, 

^ (fc)=e2 (£)/ dM^(fc)Fi(Mi) + V^(fci)F 2 (Mi)]- (92) 

In view of the explicit e 2 dependence of equation (|92p we should evaluate £k and if)±{k) in the massless limit. As 
follows from equation ([57)) 

It is readily seen from Appendix YK[ that the non-trivial kernels Fi and F2, are in fact functions of the dimensionless 
variables k and ki] note that we worked in units where H = c = 1 there. Following Damle and Sachdev^ we denote 

F t (fc, = F 2 (k,k 1 ) = $ 2 (k,h)- (94) 

The longitudinal displacement may thus be written 

e 2 f k B T^ 2 



^ =^ J \Wc) G(fc) ' (95) 



where we introduce a universal scaling function of the dimensionless variable k 



poo 

G{k) = k dh 
Jo 



■*i(Mi) + *a(Mi)- 



(e k - l) 2 v ' v ' ' (e fc i - \f 



(96) 
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We plot an appropriate moment of this distribution in Fig. [4j The heat current may be obtained from (|95[) by 
combining (|8"5|) with the usual relation (fl"7j)) 

^ efe v fc (k.U) V ,(fc). (97) 

It follows from the definitions (|17p that the corresponding thermal conductivity is given by 

_ 2c 2 f d d k 

Performing the angular integrals gives 



dTJi^W" 2 ^ (98) 



2c 2 S d 
dT (2irh) d 



dkk d+l ip{k), (99) 



where Sd = 2n d/2 /T(d/2) is the surface area of a unit hypersphere in (i-dimensions. Substituting the result (|95|) into 
((99)) and rescaling the momentum integral yields 

*—«'*■« (M Off 3 (^) r* ii+iG < E >- <'») 

Equivalcntly 

R xx =ge 2 (k B c)l%\ T {d+3 \ (101) 

where 



h he , 

are the magnetic length, and a suitable thermal wavelength respectively. The numerical coefficient is given by 

S =-J^ / dkk d+1 G(k). (103) 



d{2nY 

From equation (|101[) we see that R has dimensions Jk~ 1 m~( d ~ 2 ' > s~ x . This is consistent with the defining relations 
(|17[) . In addition, it is readily seen that our expressions for the thermal conductivity are in agreement with the scaling 
form ([8]) with dynamical exponent z = 1. To lowest order in the epsilon expansion, the integral (|103[) should be 
performed in d = 3: 

1 r°° - - 

,g = — / dkk A G{k). (104) 



37T 2 , 

This equation mirrors (3.33) of Damle and Sachdev.— The extra factor of momentum arises because we are considering 

heat transport as opposed to electrical transport. In Fig. 2] we plot the universal scaling function k A G(k). We find 
by numerical integration that 

g « 5.55. (105) 

In particular, in d = 2, we find 

Rxx = Ro (^Jp^J ' ( 106 ) 

where Rq = gk^/(Ae 2 h 3 c 4 ). That is to say, in stark contrast to the case where B = 0, the thermal conductivity does 
not diverge, but is finite^ Moreover, the dependence on e 2 is inversely related to that of the universal DC electrical 
conductivity 
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FIG. 4: Longitudinal scaling function, G(k), required for computation of the thermal conductivity in the drift regime. This 
figure was obtained by truncating the upper limit of integration at ki = 100. 



where M a ~ 1.037. When specialized to two-dimensions our result^ (jlOip may be cast in the equivalent forr 



A 



B 2 



a Q 



(108) 



where, 5 c h = 2C3fcsA T 2 , is the entropy density of a charged scalar field calculated within the epsilon expansion, and 
the transformed dimensionlcss prcfactor is given by 



A 



gK_ 

8ttC 3 2 ' 



(109) 



With g w 5.55, M a « 1.037, and C 3 = 2tt 2 /45 one obtains A « 1.19. 



X. HYDRODYNAMIC LIMIT AND INTERPOLATION 

The arguments presented above are manifestly iterative in the epsilon expansion and implicitly assume that e 2 < B. 
It was subsequently pointed out in Ref. [fusing relativistic hydrodynamic arguments directly in d — 2, that equation 
(|108p holds as an exact relation in two dimensions with A = 1, or pmhd ~ 4.66. Moreover, it was shown that all the 
physical response functions are governed by combinations of thermodynamic variables and the single parameter o~q. 
In order to make contact with the hydrodynamic results of Ref. we must study the ultra low field limit B <C e 2 . 
Even though e 2 is a small parameter, we must consider both the frequency and the field scales much smaller than 
this in order to enter the hydrodynamic regime. This is in conformity with the original zero field treatment of Damle 
and Sachde\i2£ where the hydrodynamic crossover occurs at uj ~ C(e 2 ). The magnetic field provides an additional 
cyclotron frequency and in the hydrodynamic limit this is assumed to be less than the scattering rate. As we shall 
discuss, the thermal conductivity in fact interpolates between these two closely separated limits (^4« 1.19 and A = 1) 
as the ratio B/e 2 is varied. This mirrors recent findings in strictly two-dimensional graphene where (aside from the 
important aspects of spatial dimensionality and carrier statistics) the fine structure constant, a, plays a similar role 
to eM- In the subsequent discussion we derive the exact form of the hydrodynamic result analytically, using the QBE 
and the epsilon expansion around d = 3. We will also discuss the interpolation as the magnetic field is varied. 

Returning to our original equation (|86|) . we wish to solve this coupled problem for the longitudinal and transverse 
distribution functions in more generality. It is convenient to decompose this relation in to two distinct equations 
which are respectively even and odd under reversal of the magnetic field: 



v fc .U 



e k ?^ + Q 2 B^ ± (k) 



= i4[k.u^(fc)], 



(110) 



and 



TQv t .(U x B) ifj(k) = L±[±Qk.(U x B)^ ± (k)}. 



(Ill) 
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Direct elimination will yield equations governing the longitudinal and transverse response. In order to expose this 
more clearly, we use the expressions (|A9|) and (|A26|) for the distinct longitudinal and transverse linearizations of the 
collision term: 



v fc .U 

and 



k.u|-e 2 (^)^ dfc 1 [F' 1 (fc,fc 1 ^(fc)+F^(fc,fc 1 )^(fci)]|, (112) 



T v fc .C^(fc) = ±k.C j-e 2 (£) J dh [Fi(fc, h)ip x (k) + F 2 (fc, fci)^(fci)]} , (113) 

where C = Q(TJ x B), and we have restored the factors of h and c in the collision terms. The factors of k.U and k.C 
are readily cancelled leaving coupled integral equations for the distribution functions. For pedagogical purposes it is 
convenient to discretize these equations and write them in the simpler matrix form 



where we adopt the useful shorthand 

/>oo 

M fc , fcl V i fel = / dh [Fx (k, ki)i>x(k) + F 2 (fc, k^ih)] , 



( 114 ) 



(115) 

M' kM ^ kl = I dkt [F[(k, ki)1>{k) + F' 2 (k, fci)V(fci)] • 



Rearranging the first of equations (|114|) for ?Jj± and substituting in to the second immediately yields a Fredholm 
integral equation of the second kind for the longitudinal distribution function 

i>h = i'k - (qJ-) e k M kM e kl M' kuk ^ k2 , (116) 



where 



s ekMk ^ {-§£) ■ (117) 

In the massless limit, we see that t/'oo is nothing but our previous solution obtained from equations (|92[) and (193 j) . 
It is recovered from the general integral equation (|116|1 in the limit B 3> e 2 . Employing dimensionless variables, 
k = ck/kgT, and noting the rescaling of the integration measures in (| 1 1 5|) . the solutions of equation are 
governed by the dimensionless prefactor of the second term which we may denote as 

r-i = e*(^)\ (118) 



where Ib = \/ h/QB is the magnetic length, and At = hc/ksT is the thermal wavelength. Equivalcntly, 

r~(^Pr m ) 2 , (119) 

where ~ e 2 ^T/fi is the inelastic scattering rate due to the collisions at the clean fixed point, and w* yp ~ 
QB j '{kgJ ' I 'c 2 ) is the typical cyclotron rate of a thermal carrier^ In the limit, r — > 0, the particles experience a large 
number of collisions per typical period of revolution in the magnetic field. This corresponds to the hydrodynamic 
limit studied in Refs I30ll34l In this respect, the magnetohydrodynamic parameter, r, plays a similar role to the 
Knudscn parameter in the development of Chapman-Enskog theory^ Deep in this hydrodynamic regime, B <C e 2 , 
the longitudinal distribution function satisfies the homogeneous equation 



M fejfcl £ fel M' fclifc2 V> fc2 =0, 



(120) 
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as follows directly from equation f| 1 1 6|) . In this limit, the distribution function projects on to a zero mode of the 
transverse operator M'. From our previous discussion of the transverse response (|67p it is readily seen that 

f±(k) = f {k) +Wk.U f ] + transverse (121) 



is a zero mode of this collision operator. Moreover, this form would also emerge in a relaxation time approximation 
where Af = t^c 2 and t 7 is a suitable time scale. With the parameterization (|85|) we therefore consider 

*=.V (-§£), (122 ) 

where J\f is a dimensionfull parameter to be determined. This normalization is fixed by the original inhomogeneous 
equation (|116[) . and may be found by taking moments and extrapolating to the appropriate hydrodynamic limit. 
Placing the last term of equation (|116[) on the left hand side, and inverting the matrix operation in the definition of 
ip^ ', equation (|1 16[) may be rewritten in the equivalent form 



Q 2 B 2 hc 



In order to make contact with the functional dependence (|108|) observed in the complementary regime, it is convenient 
to multiply this equation by e k and integrate over all momenta. Upon sending i\) k to the hydrodynamic form (|122p . 
and working in the massless limit where e k = ck, one then obtains the normalization condition 

ddk ,.--2,*,-U ( d fo\ 2 [ ddk ,2 ( d .fo 



where we exploit the fact that ip k is a zero mode of M'. Using equation (|15)) we see that the right hand side of this 
equation is related to the entropy density. Moreover, using equation (3.28) of Damle and Sachdev^ we see that 
the left hand side involves the DC conductivity in the absence of a magnetic field. Recalling the main steps, we 
parameterize, /± = /o ± Qk.E V'Ds(fc), and substitute into the QBE equation (fSTj) with B = 0. To linear order in the 
electric field one obtains 



c df 

-e 



2 (£)M fe , fel ^ DS (fci). (125) 



k de k 

With this identification equation (|124[) may be recast in the form 

MQ 2 B 2 c J J^ k ^ DS (k)=dT^, (126) 

where S c h is the entropy density of a charged scalar field. This may be rearranged to read 

~2Q 2 c f d d k 



MB 2 

Since the electric current is given by 



(27rft) ( 



kipBsik) 



TS ch . (127) 



J e = 2Q 2 J ^^Vfek.EVWfc), (128) 
the quantitity in square brackets is the conductivity. That is to say 

where <tq is the universal and non-trivial value of the DC electrical conductivity computed within the epsilon expansion 
in the absence of a magnetic fields From this normalization, J\f = t^c 2 , we may extract the characteristic damping 
time scale, r 7 , of the collective cyclotron mode discussed by Hartnoll et al£^ Returning to our distribution function 



(|122| . the heat current is given by 
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Oi 4.8 




FIG. 5: Evolution of the dimensionless thermal conductivity prefactor, g, as the ratio, r oc B 2 /e 4 , varies. Within the accuracy 
of our three dimensional Monte Carlo integrations (approximately 3%) the result interpolates between the value reported in our 
previous Letter,— g^ ~ 5.55, and the magnetohydrodynamic value,— go = 8-7r(27r 2 /45) 2 /(1.037) ~ 4.66. The latter corresponds 
to the exact relation, k xx (B) — TS 2 h /B 2 ag, for B < e 2 . 



The corresponding thermal conductivity reads 

2Mc 2 



dT 



d d k 



where we employ equation (|48[) again. Combining equations (|129[) and (| 1 3 1 [1 one obtains the relation 



B 2 



(131) 



(132) 



deep in the hydrodynamic limit where 6 « f 2 . We see that the functional dependence is the same in both limits, 
B 3> e 2 , and fi « e 2 , and only the dimensionless prefactor is modified by a factor close to unity. As discussed in 
Appendix [Bl numerical solution of the integral equation (|116[) yields the interpolation between these two regimes — 
see Fig. [5] It is instructive to note that the hydrodynamic result (|132p . which holds exactly in two dimensions, is 
accessible within the framework of the epsilon expansion about d = 3. The main delicate points are that we should 
work in a regime where B <C e 2 , and be careful to interpret the observables k xx , ctq, and S c h by means of their 
respective epsilon expansions. 



XI. CONCLUSIONS 



In this work we have examined the magnetothermoelectric response in the vicinity of a quantum critical point. We 
investigate the electrical and thermal transport and thermodynamics, and have presented general scaling arguments 
valid for arbitrary dimension, dynamical exponent and carrier statistics. These are supported by explicit calculations 
at the particle-hole symmetric SF-MI transitions of the Bosc-Hubbard model. The presence of a magnetic field 
strongly influences the physical response, and we demonstrate the existence of a finite thermoelectric tensor, a xy , 
and a finite thermal conductivity, R XX (B), even in the absence of impurities. We relate these observations to a 
number of different approaches, based on Lorentz invariance, the quantum Boltzmann equation (QBE), and field 
theory considerations. In accordance with recent findings of Miillcr et al for graphene,— the dimensionless prefactor 
of our thermal conductivity is a smoothly varying function of Lo c T- m . We derive an illuminating integral equation 
to describe this evolution, which exemplifies the role of zero modes in the hydrodynamic limits Our analytic and 
numerical calculations smoothly interpolate between the result presented in our previous Letter,— and the limit of 
two-dimensional rclativistic magnctohydrodynamicsi2£ The recovery of the relativistic hydrodynamic relations within 
the framework of the epsilon expansion around three dimensions is quite compelling. Although the epsilon expansion 
is well established for the calculation of critical indices relating to thermodynamic quantities, it is much less widely 
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employed in transport situations. The present body of results clearly demonstrate that we may address critical 
fluctuations in transport coefficients by using such methods. We see that the QBE approach not only has broad 
applicability but also provides a physically intuitive way to incorporate both quantum and thermal fluctuations in the 
hydrodynamic regime. Moreover, the results obtained are physically transparent and provide a platform for further 
studies. 
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APPENDIX A: COLLISION TERM AND POLYLOGARITHMS 



In order to make progress with the Boltzmann equation it is useful to linearize about an equilibrium Bose distribu- 
tion. There are two distinct linearizations depending on the context. We gather some useful formulae below. 



1. <5/±(k) = ±k.Ch(k) 



In this linearization we consider departures from equilibrium of the form <5/±(k) = ±k.C h(k) where C is a constant 
vector, and h(k) is a function of |k|. This situation arises in the longitudinal response to an electric field^ C = QE, 
and the transverse response to a temperature gradient, C = Q(U x B). We want to expand the collision term (|T5]) 
to linear order in the applied field C. We may write 



I± = - 



2ul 



where, in units where h = c = 1, 



n 



d d h 1 



(2^) d <5(k + k x - k 2 - k 3 ) (2n)5(e + e 1 -e 2 - e 3 ) 



represents the remaining phase space measure in equation (|13[) . It is readily verified that 



±ki.Ch(fci) (e fe - e ^ Ea+E3 >) n(e k )n(e 2 )n(e 3 ) 

±3k 2 .Ch(k 2 ) (e^ ~ e0O*+*O) n(e k )n(e 1 )n(e 3 ) 

±k 3 .C h(k 3 ) (e fi ^i) _ e fe ) n ( £fc ) n ( ei ) n ( ea ) + (C 2 ), 



(Al) 



(A2) 



(A3) 



where n{e) = fo(e) is the Bose distribution function, and we have used the identity 1 + n(e) = e l3£ n(e). Upon 
substituting (|A3[) into (|A1|) it is convenient to interchange k 3 <-> k 2 in the last term. This transformation preserves 
the integration measure and the linearized collision term becomes 

I±^L±=±f-_|_) Jdn [k.C/i(fc)T 1 (fc 1 ,A :2 ,A :3 )+k 1 .C^(fc 1 )T 2 (fc,A :2 ,A :3 )+k 2 .C/ l (A :2 )T 3 (fc,A :i ,A :3 )] (A4) 
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where 

Tifai.fc.fca) = 3(e^ 2+£3 )- e ^ 1 )n( £l )n( e2 )n(e 3 ), (A5) 

T 2 (fc, fe, fc 3 ) = (e** - e ^ 2+ ^) n(e fe )n( £2 )n( £3 ), (A6) 

T 3 (Mi,fc3) = 2(e^ 3 -e^ £ * +£l ))n(e fc )n(ei)n(e 3 ). (A7) 

As discussed by Damle and Sachdev^ to leading order in the epsilon expansion one may evaluate the necessary 
integrals directly in d = 3. In addition, one may consider the massless limit where = k, and we have set c = 1. 
The angular inte gral s and one of the radial integrals may be carried out explicitly, by means of formulas (CI), (C3) 
and (C5) of Ref. l20t in their notations cPki = kfdkidtli. Interchanging ki <-> k 3 in the first term, and ki <-> k 2 in 
the last, this procedure yields 

/ 2u 2 \ kC r°° ( h(k) 

L± = ± (^-^J (4^)3 J dhdk 2 l-^T 1 (k 1 + k2-k,k 2 ,h)l 1 {k,k 1 ,k 2 ) 

} [T 2 (fc, k 2 ,k + h- ha) l 2 (k, kt,k 2 ) - T 3 (k, k 2 ,k + k 2 - h)l 3 (k, k x , ha)]] , (A8) 



3fc 4 

where Ii, I 2 and I 3 are (domain dependent) polynomials given in (C2), (C4) and (C6) of Ref. [2(| In this way one 
may write the linearized collision term in the compact form2£ 

L± = ±k.c|-e 2 ^ dki [h{k)Fi(k,ki)+h(ki)F 2 (k,ki)]^ , (A9) 

where we use the fact that uq = (24/5)7r 2 e, and we define 

Fi(ft,*i) = J dk 2 Ti{hx + k 2 -hMM)HKki,k 2 ), (A10) 

and 

F 2 (k,k 1 ) = ~y^4 / dk 2 [T 2 (k,k 2 ,k + ki-k 2 )I 2 (k,ki,k 2 ) 

-T 3 (k, k 2 ,k + k 2 - fc 1 )I 3 (fc, k u k 2 )\ . (All) 
We again emphasize that the functions Ti, T 2 and T 3 are evaluated in the massless limit. More explicitly 

Fi(Mi) = — £5^/ dfen(*2)[l + n(*i + *2-*)]Ii(fc,Ai,fe)- (A12) 

Likewise, if we denote 



then 



F 2 (fc, h) = Ff (fc, *i) + F§(fc, fci) (A13) 
Fa(fc.*i) = ^ [1 J"?f fc ? ] / dfca»(*aK* + fci-*a)Ia(fc,Ai,*a), (AH) 

47T 77 fA-l f°° 

F b 2 (k,h) = -— ^> J d*an(*a)[l + n(fc + *a-Ai)]l3(fc,*i,fca). (A15) 

These expressions are in conformity with equations (C7) and (3.28) of Ref. [2(J As noted by Damle and Sachdev. the 
integrals (|A12[) . (j A14|) and (|A15[) may be evaluated exactly using polylogarithm functions. Although the method was 
carefully explained, the explicit form of these kernels was not stated in their original works . 20 ' 21 The expressions are 
quite lengthy, and are rather tedious to derive. Nonetheless, they are valuable for numerical work. We may write 

Fi(fc, fci) = l "^ 1 ,^,: kl) [6(fc - h)tx 2 (k, h) - 6(fci - k)(x 2 (h, k)) , (A16) 
25 k n(K) 
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together with 



and 



F a 2 (k,h 



2tt [l + n(fc)]n(fc + fci) fc , 
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F 6 U- J- ^ - 4?r - fc) 



fc 4 n(fci) 



[e(jfc-*i)L|(Ar,jfci)-e(ifci-*r)LS(ibi,fc)] , 



where 0(fc — fci) is the step function. In writing these kernels we have introduced 

LS(fc,fci) = 24A4(fc,*i) + 12[fcj7 3 (A! ) fci) + An73(fci,*)]-6fcfciA3(* J *!i)i 



where 



Further. 



where 



A±(x,y) ee /T" [Li n (e-^) + Li n (e-^)±Li„(e-^ a:+I '))±Li„(l) 
r; n (rE,y) ee /T n [Li„(e"^) - Li„(e"^) - Li^e""^) + Li„(l) 

L^(fc, fcx) = -3 [4^4 + 2(fc - fci)/ti 3 - fcfci/12 + 4fci^ 3 + 2fcfc x i^] , 
/i n (x,2/) == /T" [Li„(l) + Li„(e"^) - Li„( e -^) - Li^e"^-^) 



and 



i/ n (ar,y) ee /T" [Li„(e-^) - Li n (e-*0] 
Here Li p (z) is the polylogarithm with series expansion 



u p( z ) = 



z" 
■nP' 



(A17) 



(A18) 



(A19) 



(A20) 



(A21) 



Note that we have used the fact that \j 2 {k,k\) is symmetric in order to eliminate the step functions from Frj^fc, fci). 
In writing these expressions we employ polylogarithms whose arguments lie within the unit disc. Although tedious 
to check analytically, the equality of these functions and the integral representations (|A12[) . (|A14|) and (|A15[) is 
readily verified numerically. It is worth noting that the kernels Fi(fc, fci) and F 2 (k, fci) possess singularities when their 
arguments coincide. 



2. <5/±(k) = k.Cfc(fc) 

In this linearization we consider departures from equilibrium of the form <5/±(k) = k.C h{k) where C is an arbitrary 
vector, and h(k) is a function of |k|. This charge independent situation arises in the longitudinal response to a 
temperature gradient, C = U ee (— VT)/T, and in the transverse response to an electric field, C = E x B. In this 
case 

^ ut - J*£ = 3k.C h(k) f e 0(«*+«) - e^ 1 ) n(e 1 )n(e 2 )n(e 3 ) 

+ 3k!.C/i(fci) ( e ^+ £ 3) _ e /fc») n(e k )n(e 2 )n(e 3 ) 
+ 3k 2 .C h(k 2 ) (e^ - e « e *+ ei >) n(e fc )n(ei)n(e 3 ) 

+ 3k 3 .C h(k 3 ) {e^ - e ^*+ 6 ^ n(e k )n( £l )n(e 2 ) + 0(C 2 ). (A22) 
Substituting this into the collision term, and again making the interchange k 3 <-> h 2 in the last term, we find 

I±^L' ± =(-^\ J dn [k.C / 1 (fc)Ti(fc l! fc 2 ,A ;3 ) + k 1 .C h(k 1 )T' 2 (k,k 2 ,k 3 ) + k 2 .C h(k 2 )T 3 (k,h, fc 3 )], (A23) 
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where 



Ti(fci, k 2l ka) = 3 (e^+^ - e^ 1 ) n(e 1 )n(e 2 )n(e 3 ) 
T' 2 (k,k 2 ,k 3 ) = s(e^+^ -e^n(e k )n(e 2 )n(e 3 ) 
T^(fc,fc 1; fc 3 ) = 6(e fe -e^ +£l ))n( £fc )n( £l )n( e3 ). 
In particular, it is readily seen that Tj = Ti, T' 2 = — 3T2 and T 3 = 3T 3 . This procedure therefore yields 



(A24) 



14 = I - 



2u2\ k.C f°° „ „ (h(k) „ 



+ k 2 - k, k 2 , ki) li(k, ki,k 2 ) 



t T 2( fc ' fc 2, k + ki — k 2 )h(k, k u k 2 ) + T 3 (fc, k 2 , k + k 2 - fci)I 3 (fc, k u k 2 )] 



That is to say, the linearized collision terms may now be written 



L' ± = k.C < 



dk x [F[{k, ki) h(k) + F' 2 (fc, fa) h(ki)] } , 



where 



Fi(k, fa) = F 1 (k, fa), F' 2 (k, fa) = 3 (F b 2 (k, fa) - F"(fc, fa)) , 
and Fi(k,fa), F^(k,fa) and F\{k, fa) are given by equations (|A16|) . (|A17|) , and (|A18|) respectively. 



(A25) 



(A26) 



(A27) 



APPENDIX B: LONGITUDINAL CROSSOVER EQUATION 

As discussed in section |Xl the thermal conductivity is a function of the dimcnsionlcss parameter, r oc _B 2 /e 4 , 
which controls the ratio of the typical cyclotron frequency for a thermal carrier to the inelastic scattering rate. The 
longitudinal distribution function satisfies the integral equation ()116[) 



-(- 



\QBhc 



e k M kM e kl M' kiM ip k2 , 



(Bl) 



where ip™ , M and M' are given by equations (|117[) and (| 1 1 5[) respectively. For both numerical and analytic purposes 
it is convenient to recast this equation in terms of the dimensionless momenta, k = ck/ksT. This rescaling modifies 
the coefficient of ip^ given in equation (|117[) . and it is convenient to introduce 



(k B T 



h \QBc 



vf(fc), 



(B2) 



and similarly for Adopting this rescaling the integral equation (|116[) may be recast in the dimensionless form 



*(fc)=* 00 (fc)-r- 



dfc 1 [Q 1 (fc,fc 1 )*(fc) + Q 2 (fc,fc 1 )*(fc 1 )], 



where 



and 



(QB\ 



2 ' hc^ 4 



e 4 V h J \ k BT 



*oo(fc) = k I dki 
Jo 



ke 



(e k - 1) 



;<f>l(k,kl) + $2{k,kl) 



he*" 1 
(e*i - l) 5 



(B3) 



(B4) 



(B5) 
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coincides with our previous distribution function (|96[) . The non-trivial kernels, Qi(fc, ki) and Q2(fc, k\), are given by 
the integrals 



Qi(fc, fci) 



(B6) 



and 



Q a (ft,*l)= / *2 [fc$l(fc,fc 2 )fc$2(fc^l)+fc $ 2(fc,fcl)fcl$i(fcl,fc 2 ) + fc$2(fc,fc2)fc2*2(fc2,fcl)] , (B7) 

Jo 

where F^fc, k{) = $i(fc, fci), F 2 (fc, fei) = $2(^1 and their primed counterparts, are the dimensionless kernels given 
in Appendix [A] Note that in deriving the results (|B5j) , (|B6[) and (|B7|) , we have also taken the massless limit in 
accordance with the epsilon expansion. 

It is evident from the longitudinal crossover equation (|B3|) . that the distribution function interpolates between 
^oa(k), as r — > 00, and a zero mode of the integral operator in the hydrodynamic limit, r — > 0. In order to see this 
more formally it is natural to consider an expansion of the form 



(fc) = *o(£) + + r 2 * 2 (fc) 



(B8) 



Substituting this expansion into equation (|B3|) and equating coefficients at order r 1 one finds that for a non- vanishing 
^o(fc) to be present it must be an exact zero mode of the integral operator: 



dh[Qi(k, fci)^o(fc) + Q 2 (fc, h^oih)} = 0. 



As discussed in section [Xl the function 



Wo(kJ oc — — oc 



de k (e k - If 



= n (k), 



(B9) 



(BIO) 



satisfies this homogeneous condition. Although it will not concern us here, equating coefficients at higher order in 
r leads to a recursive hierarchy of integral equations for the functions ^i(k). These are similar in spirit (if not in 
details) to those encountered in the Knudsen expansion of the Boltzmann equation; see for example §6.2 of the book 
by Harris^ Equation (|B9[) is particularly important from a numerical perspective, since any regularization which 
lifts this zero mode property (e.g. through truncation of the integration limits or through rounding errors) may 
potentially yield a solution starting at order r, instead of order r°, as evident from equation (|B8jk on dimensional 
grounds this would yield a thermal conductivity of the Wiedemann-Franz form. In order to recover the exact non- 
vanishing hydrodynamic limit (| 1 32[) . it is essential that this zero mode feature is properly implemented. To this end, 
let us parameterize our solutions to the integral equation in the form 



= 72o(fc)$(fc), 



(Bll) 



where $(fc) is the solution to be determined. Assisted by the exponential decay of the prefactor, we convert the 
integral equation (|B3[) into an approximate matrix equation by expanding $(fc) in a basis of Nf, basis functions: 



N h -1 



(B12) 



n=0 



Substituting this decomposition into equation (|B3|) . multiplying by k 2 f m (k) and integrating (where we incorporate 
an extra factor of k 2 for convergence purposes) one obtains the matrix equation 



where, A r 



and 



A 1 



L Aj^ n , with explicit matrix elements 



A 1 = 

ran 



dk k 2 Tlo(k)f m {k)f n {k), 



A 



11 _ 



dkdhdhk 2 f m (k) n (^f n (m[ A \k,h,h) + QT\k,k 1 ,k 2 )n (h)f n (h) 



(A), 



(B13) 



(B14) 



(B15) 
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where Q^ A \k, k\, ^2) denote the integrands (or arguments) of the integral representations (|B6[) and (|B7[) . In addition 

/•OO 

B m = / MlP/mW [%(^l(Ml) + $2(Ml)Wfe)] • (B16) 

Jo 

The dimensionless thermal conductivity parameter is defined as 



9=^2 I 



and is thus approximated as 



1 N J'~, 1 r°° 

n ^c„ dkk 4 n (k)f n (k). (bis) 



3tt 2 

n=0 



Evaluating the matrix elements numerically (whilst implementing the zero mode condition Aj^ = exactly) we may 
solve the linear system of equations for the coefficients c„. Using a basis of iVj, = 7 monomials (1, k, . . . , k 6 ) we plot 
g as a function of r in Fig. [5] It is readily seen that within the numerical accuracy of our Monte Carlo integrations 
(approximately 3%) this dimensionless coefficient interpolates between the result reported in our Letter^ and the 
hydrodynamic result^ as the parameter r oc B 2 /e 4 is varied. 
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In general it is necessary to distinguish between the microscopic currents (| 15p and (|16p . and the transport currents (|17l) . 
due to the presence of magnetization currents . 41 ' 56 This issue is circumvented in conventional Boltzmann approaches to 
transport in magnetic fields which do not include magnetization corrections to the distribution function. 
Note that at weak fields, the B 2 contributions to the free energy are eliminated in these works by a renormalization of the 
electric charge associated with the Maxwell term. 

We are extremely grateful for illuminating conversations with David Huse on this point. 



